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ABSTRACT 

Context. To understand the origin of stellar activity in pre-main sequence Herbig Ae/Be stars and to get a deeper insight in the interior 
of these enigmatic stars, the pulsational instability strip of Palla and Marconi is investigated. In this article we present a first discovery 
of non radial pulsations in the Herbig Ae spectroscopic binary star RS Cha. 

Aims. The goal of the present work is to detect for the first time directly by spectrographic means non-radial pulsations in a Herbig 
Ae star and to identify the largest amplitude pulsation modes. 

Methods. The spectroscopic binary Herbig Ae star RS Cha was monitored in quasi-continuous observations during 14 observing 
nights (Jan 2006) at the Im Mt John (New Zealand) telescope with the Hercules high-resolution echelle spectrograph. The cumulated 
exposure time on the star was 44 hrs, corresponding to 255 individual high-resolution echelle spectra with R = 45000. Least square 
deconvolved spectra (LSD) were obtained for each spectrum representing the effective photospheric absorption profile modified 
by pulsations. Difference spectra were calculated by subtracting rotationally broadened artificial profiles; these residual spectra were 
analysed and non-radial pulsations were detected. A subsequent analysis with two complementary methods, namely Fourier Parameter 
Fit (FPF) and Fourier 2D (F2D) has been performed and first constraints on the pulsation modes have been derived. 
Results. For the very first time we discovered by direct observational means using high resolution echelle spectroscopy non radial 
oscillations in a Herbig Ae star. In fact, both components of the spectroscopic binary are Herbig Ae stars and both show NRPs. 
The FPF method identified 2 modes for the primary component with (degree azimuthal order m) couples ordered by decreasing 
probability: /, = 21.11 d"' with (Cm) = (11,11), (11,9) or (10,6) and = 30.38 d"' with (Cm) = (10,6) or (9,5). The F2D analysis 
indicates for fi a degree f = 8- 10. For the secondary component, the FPF method identified 3 modes with (C,m) ordered by decreasing 
probability: /, = 12.81 d"' with (t,m) = (2,1) or (2,2), /z* = 19.11 d"' with (Cm) = (13,5) or (10,5) and /a = 24.56 d^' with {i,m) 
= (6,3) or (6,5). The F2D analysis indicates for fi a degree f = 2 or 3, but proposes a contradictory identification of /2 as a radial 
pulsation (f = 0). 

Key words, stars: pre-main sequence - stars: oscillations - stars: individual: RS Cha - stars: binaries: spectroscopic 



1. Introduction 

Asteroseismology represents a modern tool for studying 
the stellar interiors of the enigmatic group of pre-main se- 
quence (PMS) stars of intermediate mass (~ 2-8 M©), so 
called Herbig Ae/Be stars ( [Herbig, 1960[ [Strom etal., 1972[ 
[Finkenzeller & Mundt, 1984[ iFinkenzeller & Jankovics, 1984|l' 
Since their first systematic classification in 1960 (Herbig, 1960 1 
this group of young stars has been extensively stud- 
ied, but one of the major questions remains unan- 
swered: how can the intense stellar activity and strong 
stellar winds (|Praderie et al., 1982[ [Catala et al., 1986[ 



Catala & Kunasz, 1987[ [Bohm & Catala, 1995| and 
Bohm et al., 1996 1, as well as the many highly variable emission 



lines observed in the spectra of these stars be explained? 

Magnetism is frequently invoked as being responsible for 
active stellar phenomena. The position of the Herbig stars in 
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* Based on observations collected at the Im McLellan telescope at 
Mt John, NZ 



the HR diagram indicates that they are in the radiative phase 
of their contraction towards the main sequence ( |Iben, 1965) 
[Gilliland, 1986| l, and should in principle not possess any outer 
convective zone. In its absence, the classical magnetic dynamo 
mechanism could not be responsible for the observed active phe- 
nomena. Moreover, recent spectropolarimetric observations of 
Herbig stars indicate the presence of significant large-scale mag- 
netic fields only in a small fraction of them ( [Wade et al., 2005[ 
[Catala et al., 2007] and [Alecian et al., 2008| l. This fraction, close 
to 10%, is consistent with a magnetic field distribution similar to 
that of the main sequence stars. In that scenario, magnetic Herbig 
Ae stars would be the progenitors of main sequence Ap/Bp stars. 
These recent results are in agreement with the primordial fossil 
field hypothesis. Unless magnetic fields in Herbig Ae/Be stars 
are of very complex nature and can therefore not be detected 
by spectropolarimetric means, in opposition to clearly organized 
dipolar fields, the pure magnetic origin of activity in these stars 
is more than doubtful. 

For many years non-classically generated magnetic fields 
have been used as a potential explanation for the observed ac- 
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tive phenomena. Even if the presence of complex fields can not 
be ruled out it is of major importance to investigate other pos- 
sible external or internal origins of this tremendous amount of 
dissipated energy, as witnessed by chromospheres and coronae, 
but also variable spectral lines, winds and bipolar jets. 

The only way of studying in detail the internal stellar 
structure is the analysis and the modelling of stellar pulsa- 
tions, if observed. As an example, PMS stars gain their en- 
ergy from gravitational contraction and therefore differ signi- 
cantly from post-main-sequence stars having already processed 
nuclear material - stellar pulsations are sensitive to these dif- 
ferences of internal stellar structure (see e.g. |Suran et al., 2001| l. 
As of today, the internal stellar structure of PMS stars is not 
yet well constrained. Since few years the existence of pul- 



Table 1. Parameters 
lAlecian et al., 2005[ 



sating intermediate mass PMS stars is known ( Breger, 1972 
Kurtz & Marang, 1995] [Donati et al., 1997| l. This observa tional 



result motivated Marconi & Palla ( Marconi & Palla, 1998) 1 to in- 
vestigate the pulsation characteristics of HR 5999 theoretically, 
which enabled them to predict the existence of a pre-main- 
sequence instability strip, which is being crossed by most of 
the intermediate mass PMS objects for a significant fraction of 
their evolution to the main sequence. This strip covers approxi- 
mately the same area in the HR diagram as the S Scuti variables. 
Zwintz ( jZwintz, 2008| l compared, based on photometry, the ob- 
servational instability regions for pulsating pre-main sequence 
and classical 6 Scuti stars and concluded that the hot and cool 
boundaries of both HR diagram instability regions seem to co- 
incide. This preliminary result deserves further study by aim of 
full asteroseismological approach based on spectroscopy. 

As of today, more than 30 intermediate-mass PMS stars 
have revealed to be pulsating at time-scales typical of 6 Scuti 



stars (see e.g. Kurtz & Marang, 1995 


IKurtz & Catala, 2001 


Donati et al. , 1 997| Bohm et al. , 2004 


Marconi et al., 2002 


Ripepi & Marconi, 2003 |Zwintz & Weiss, 20031 ICatala, 2003 



and references therein). 

RS Chamaeleontis is a bright spectroscopic eclipsing binary 
star. Both components are Herbig Ae PMS stars of similar mass 
(close to 1.9 Mq). Recently the age of RS Cha has been de- 
termined to Myr (Luhman & Steeghs, 2004 1, which verifies 
it's PMS nature. [Andersen, 1975 [ already reported small ampli- 
tude radial velocity variations on top of the binary radial ve- 
locity curve for both components of RS Cha, suggesting the 
possible presence of stellar pulsations. Photometric observations 



by Mclnally & Austin, 1977 revealed short-term variations in at 
least one of the two components, possibly linked to stellar pulsa- 
tions. Very recently, Alecian et al. ( [Alecian et al., 2005] l reported 
radial velocity variations in the residual velocity frame (cleaned 
for orbital velocity) with amplitudes up to a few km s ' and pe- 
riods of the order of Ih, indicative of 6 Scuti type pulsations. 

The aim of our study of the two components of RS Cha is to 
provide a first set of asteroseismic constraints for forthcoming 
non-radial pulsation models by determining unambiguously a 
higher number of periodicities and identifying, in a second step, 
the corresponding pulsation modes with their respective degree 
{ and azimuthal number m. 

To achieve this goal, we decided to perform high resolution 
spectroscopic observations on a large time basis and with opti- 
mized time coverage. 

Section |2] reviews previous related work. Sect. [3] describes 
the observations and data reduction. Sect. |4] summarizes results 
of the orbit determination. Sect. |5] reveals the detection of non- 
radial pulsations in both components of RS Cha, Sect. |6] and |7] 
present frequency analysis and moment identification in the pri- 



,Clausen & Nordstrom, 1980[ 



of R S Cha. Referenc es : [ 1 ] 
[2] [Ribas et al., 2000| [3] 



Parameter 



Primary Secondary References 



Teff [K] 

log(L/Lo) 
log(^) [cm s-^] 
V sin I [km s"'] 

Porb [d] 

i [deg] 
[Fe/H] 



1.89±0.01 
2.15±0.06 
7638±76 
1.15±0.09 
4.05+0.06 
64+6 



1.87+0.01 
1.87+0.01 
7228+72 
1.13+0.09 
3.96+0.06 
70+6 



1.67 
83.4+0.3 
0.17+0.01 



[1] 
[1] 
[2] 

g = MGIR^" 
[1] 
[1] 
[3] 
[1] 



Table 2. Log of the observations at Mt John Observatory, 
NZ, in Jan 2006. (1) and (2) Julian date (mean observation) 
(2,450,000+); (3) Number of high resolution RS Cha spectra; 
(4) typical range of S/N (pixel"') at 550 nm (centre of V band) 



Date 


JDflist 


JDla„ 


N 

^ ' spec 


S/NRange 


(1) 


(2) 


(3) 


(4) 


(5) 


Jan 09 


3745.0101 




1 


120 


Jan 10 


3745.9544 


3745.9733 


2 


80-100 


Jan 12 


3747.9537 


3748.1779 


28 


70-90 


Jan 13 


3749.0013 


3749.1930 


23 


100-120 


Jan 14 


3749.8911 


3749.9911 


12 


70-90 


Jan 15 


3750.8921 


3751.1919 


31 


150 


Jan 16 


3751.8958 


3752.1970 


29 


90-120 


Jan 19 


3754.9015 


3755.2006 


33 


100-130 


Jan 20 


3755.9245 


3756.1993 


25 


90-120 


Jan 21 


3756.9068 


3757.0386 


14 


60-150 


Jan 22 


3757.9004 


3758.1621 


32 


120-170 



mary and secondary component, respectively. A discussion is 
proposed and a conclusion is drawn in Section |8] 



2. Previous related work 

The pre-main sequence spectroscopic eclipsing binary RS 
Cha has been studied extensively throughout the last years. 
Thanks to its eclipsing nature and the known inclination an- 
gle the system has fully been calibrated ( Alecian et al., 2005] 



[Alecian et al., 2007a| [Alecian et al., 2007b[ ). Table [T] summa 
rizes the main results. 



3. Observations 

The analysis presented in this paper is based on a 14 nights 
observing run in January 2006 at the Im Mt John telescope 
equipped with the Hercules echelle spectrograph. We obtained 
quasi-continuous single-site observations of the target star dur- 
ing these 2 weeks and obtained a total of 255 individual stel- 
lar echelle spectra, each spectrum having an individual exposure 
time of lOmin. The star was observed in high resolution spec- 
troscopy at R ^ 45000 and covering the wavelength area from 
457 to 704 nm, spread over 44 orders. The detector was a Ikxlk 
Site CCD. The highest S/N (pixel"') values we obtained reached 
210 on Jan 16th, corresponding to almost 300 per resolved ele- 
ment (2 pixels); typical values of S/N (pixel"') ranged around 
80-150 in this run. Table |2] summarizes the log of the observa- 
tions. 

The general observing strategy was to obtain as many 10 
minute observations of the target star during the night as pos- 
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sible. Before the beginning and after the end of the night we 
obtained complete sets of caHbrations: in general 6 bias images, 
6 tungsten light flat field spectra, and, the closest possible to the 
first stellar spectrum, 2 ThAr arc spectra (the calibration set at 
the end of the night being scheduled in the opposite order). In 
order to guarantee the possibility to monitor directly small spec- 
tral shifts we decided to interleave regularly ThAr arc spectra 
during the night (every two hrs). 

Most of the data reduction was carried out following 
standard reduction procedures using the "ESPRIT" spectro- 
scopic data reduction package (Don ati et al., 1997 |. This pack- 
age also makes use of the "optimal extraction algorithm" (Horne, 
|Horne, 1986[ ). Heliocentric velocity and Julian date correction 
was performed. The intrinsic wavelength calibration accuracy 
achieved with the "ESPRIT" 2D-polynomial fit procedure is 
better than 0.22 pm mean rms (i.e. 2.2 mA, corresponding to 
120 ms ' at 5500 A) for the Hercules data set. 

The best way to correct the 2D-wavelength polynomial for 
evolution during the night due to the intrinsic spectral instability 
of the instrument revealed to be the computation of least square 
deconvolved LSD profiles ( |Donati et al., 1997] ) of the more than 
100 telluric water vapour lines contained in each stellar spectrum 
and to cross coiTelate them with respect to the first spectrum 
of the night. The precision of the above correction is around 
50 ms ' (see Appendix A in |Bohm & Hirth, I997| ). AU spectra 
were then corrected with their individual time-dependent correc- 
tion to the initial wavelength calibration based on the ThAr arc 
lamp spectrum acquired at the beginning of each night. More 
details about this procedure can be read in |Bohm et al., 2004] 

The next step of the data reduction was to calculate for all 
255 stellar spectra photospheric LSD-profiles, using a mask cor- 
responding to the spectral type A7, but taking care to elimi- 
nate from the line list all transitions showing non-photospheric 
behaviour (Hydrogen Balmer lines. He I D3, many Fell lines, 
Na I D,...). The equivalent photospheric profiles contain eventu- 
ally multiplex information of 1930 individual photospheric lines, 
most of them being very weak lines. Typically, a spectrum of RS 
Cha with S/N = 150 (pixeL') at 5500 A yields a typical S/N 
value of 1400 per velocity step (or pixel) in the resulting LSD 
profile. The accuracy of the radial velocity determination on an 
individual LSD profile scales, to the first order, with the inverse 
of the square root of the weighted multiplex gain and is there- 
fore estimated to be as low as 12ms"'; however, the limiting 
factor in the precision of the overall radial velocity correction 
is clearly the spectra-to-spectra cross correlation of telluric line 
LSD profiles. We therefore estimate the finally achieved preci- 
sion in radial velocity to be around 50 ms ' for each individual 
stellar LSD profile of this data set, which is largely enough for 
our spectroscopic study. 

Fig-fflshows a typical LSD profile of RS Cha with both pri- 
mary and secondary component well separated in radial velocity. 



4. Orbit determination with circularized binary 
approximation 

The spectroscopic binary LSD profile was then fitted with a dou- 
ble rotational profile (see e.g. [Gray, 2005 1 with an IDE proce- 
dure and yielded for all cases except the fully merged profiles 
precise heliocentric radial velocities of the two components. In 
a second step, the binary orbit was then fitted with a circularized 
orbit (i.e. sinusoidal radial velocity curves) on the remaining 180 
"binary"-LSD profiles and can be seen in Fig.|2] Orbital parame- 



Simultaneous fit of primary arid secondary profile 




elocity (km/s) 



Fig. 1. This figure shows a binary spectrum from Jan. 12* 2006 
(A HJD - 3745.9575) (thick line) superimposed with a simulta- 
neous rotational fit of both components. The little bumps in the 
core of the line profiles are indicative of non-radial pulsations. 



ters are listed in Table[3] The less massive secondary component 
was identified based on its slightly larger orbital amplitude. 



Table 3. Orbital parameters of both components of RS Cha. 
AHJD = HJD - 2450000.0. 



Parameter 


primary component 


secondary component 




3750.5044 


3749.6695 


Period (days) 


1.66965± 0.00003 


id 


Amplitude (kms"') 


136.4159 


137.9077 



It was interesting to determine the orbital period and its 
phase shift with respect to the work published previously. Based 
on the method described in Alecian et al. ( |Alecian et al., 2005) 1 
we calculated the observed minus computed (O-C) timings of 
the first conjunction and added our new point to the Figure [3] 
showing it's evolution over a time span of 30 years. We can ob- 
serve that the curve seems to decline slightly in gradient, indica- 
tive of either the presence of a third body in the RS Cha system 
or of the signature of changes in the orbital period of the order of 
AP/P = (7.5 + 0.7) 10"^. Further observations in the next years 
would be required to better characterize these orbital variations 
and to understand their origin. 



5. Detection of non-radial pulsations in both 
components 

Finally, we extracted the two binary components of each LSD 
profile along its individual orbital radial velocity in order to ob- 
tain two files with centred dynamical spectra revealing only the 
intrinsic stellar profile variations. Figs.|4]and|5]reveal the devia- 
tions from the mean (rotational) profile for three different nights. 
It can be clearly seen that bumps move through the line profile 
in a quite complex manner Fig. |4] (primary) reveals a more diag- 
onal trend of the bump movement (bottom-left to up-right) than 
|5] (secondary) . As will be later shown, these features are due to 
several high-degree non-radial pulsation modes. This represents 
a first direct detection of non-radial pulsations in Herbig Ae stars 
by spectroscopic means. 
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Fig. 2. Orbital movement of the spectroscopic binary components of RS Cha. AHJD - HJD - 2450000.0. The continuous and dotted 
hnes show the adjusted circular orbits of the primary and secondary component, respectively. The dots represent the observed radial 
velocity values of both components. Dot sizes do not represent error bars and have been selected only for better visibility of the 
figure. Typical error bars are of the order of 50 ms ' and can not been represented here, due to the large scale of the orbital amplitude. 
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Table 4. Finally selected data subset used for the frequency 
analysis. (1) Night and (2) Julian date of the first obser- 
vation (2,450,000-f); (3) Julian date of the last observation 
(2,450,000-H); (4) Number of spectra; (5) Duration (hrs) 



45000 
HJD - 2400000 



Date 


JDflist 


JDu,„ 


N 

^ ' spec 


Duration 


(1) 


(2) 


(3) 


(4) 


(5) 


Jan 12 


3747.9537 


3748.1607 


26 


4.97 


Jan 13 


3749.0013 


3749.1099 


13 


2.60 


Jan 14 


3749.8911 


3749.9411 


6 


1.20 


Jan 15 


3751.0417 


3751.1838 


16 


3.41 


Jan 16 


3751.8958 


3752.1861 


28 


6.97 


Jan 18 


3753.9037 


3754.1257 


16 


5.33 


Jan 19 


3754.9015 


3754.9259 


8 


1.32 


Jan 20 


3756.0609 


3756.1912 


16 


3.13 


Jan 22 


3757.9004 


3758.1621 


32 


6.28 


Total 






161 


35.2 



Fig. 3. Figure from [Alecian et al., 20031 One new full circled 
point was added, corresponding to this new 2006 data set. 

The colour coding in Figs. |4] and |5] represents a dynamical 
range of about 0.5% of the continuum. The first reference spectra 
of each night can be seen in Table[2] 

It must be noted at this stage that low frequency remnants of 
unproper orbital fitting subsist which can not easily be accounted 
for with a circularized binary orbit fit, revealing perhaps the pres- 
ence of an additional yet undetected multiple component. 

6. Pulsational analysis of the primary component 

A detailed search for frequencies and the subsequent identifi- 
cation of the pulsational modes is described in this section. At 
this stage, the spectra have been corrected for the radial veloc- 
ity variations caused by the binary orbit of RS Cha. We selected 
only spectra for the frequency analysis and mode identification 
that do not show any blending of the LSD profiles of the two 
components, i.e., the situation when the difference in orbital ve- 



locities was larger than approximately the sum of the two pro- 
jected rotational velocities, Av > vi sin ; + V2 sin /. The useful 
Doppler velocity range for the analysis was thus between -55 
and +55 km s ' which excluded the continuum. A set of 161 
spectra covering 35 hours was finally selected for the analysis. 
See Table|4]for a log of the observations. 

6.1. Frequency analysis of the primary component 

We have used the software package FAMIAS ( |Zima, 2008| l for 
the search for periodicities in the line profile and for the iden- 
tification of the pulsation modes using the Fourier parameter 
fit (FPF) method ( |Zima, 2006| l. The search for periodicities was 
carried out by using Discrete Fourier Transformation (DFT) and 
multi-periodic non-linear least-squares fitting (LSF), a method 
which fits a sum of sinusoidals to the variations. Both methods 
assume sinusoidal variations, which in general is a good approx- 
imation for pulsation modes having low amplitudes. 
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Fig. 4. Line profile variations due to the oscillations of RS Cha. Fig. 5. Line profile variations due to the oscillations of RS Cha. 

The deviations from the mean intensity are displayed for three The deviations from the mean intensity are displayed for three 

different nights (Jan 12*, 16* and 22"'' 2006). Primary compo- different nights (Jan 12*, 16* and 22"'* 2006). Secondary com- 

nent of RS Cha. ponent of RS Cha. 



In order to search for pulsation frequencies we applied the 
following iterative procedure: in the Fourier spectrum, the fre- 
quency peak having highest amplitude is selected and the data 
pre-whitened with a least-squares fit that included this peak and 
all previously found frequencies. If several peaks have similar 
amplitudes, e.g., due to aliasing, the peak which results in the 
lowest residuals is selected. The search is stopped if none of the 



remaining peak reaches a S/N of above 4 ( Breger et al., I993] l, 
corresponding to a 99.9% confidence level. 

Since the dispersion range that we selected for analysis did 
not extend towards the continuum, we were not able to prop- 
erly compute line moments. We therefore computed the centre of 
gravity (CoG) of the line profile, which is numerically similar to 



T. Bohm et al.: Discovery of non-radial pulsations in the spectroscopic binary Herbig Ae star RS Cha 



6 

the first moment, to study its temporal variations. Furthermore, 
we examined the pixel-intensity-variations (PIV) across the line 
profile for the frequency search. We studied the PIV by com- 
puting the mean Fourier spectrum across the line profile. The 
significance of a frequency peak was determined by computing 
its S/N in the wavelength bin where it has largest amplitude. The 
two methodologies are complementary in the way that the CoG 
is more sensitive to low-degree pulsation modes {£ < 4), whereas 
the PIV are also sensitive to modes of higher degree. 

The Fourier spectrum of the CoG time series showed many 
peaks in the low-frequency domain between and 6 d ' (see 
Figure |6l). These peaks are associated to nightly zero point ve- 
locity shifts of the line profile and are likely caused by the uncer- 
tainty of the binary orbit fitting. No significant peaks were found 
above 6 d"' . 

The frequency analysis of the PIV also revealed several 
peaks in the low frequency domain (see Figure|7]). Here, they are 
clearly associated with the uncertainties of the orbit fitting. We 
found significant frequencies at 0.599 and at 1.197 d"' , which 
corresponds to the orbital frequency and twice thereof After 
pre-whitening with these terms, two other significant frequencies 
could be detected in the data: /i - 21.11 and /2 - 30.38 d""' . 
After further pre-whitening, no additional frequencies having a 
S/N above 4 could be found. Table 16.1] lists the frequencies de- 
tected in the primary component. 



Table 5. Frequencies detected in the pixel-by-pixel intensity 
variations of the primary component of RS Cha. The value Apiv 
is the integral of the pixel-amplitude across the line profile be- 
tween -55 and 55 km s"' . For each frequency, the S/N is com- 
puted for the wavelength bin where its amplitude is highest. 



ID 


freq. 


Apiv 


S/N 




d-' 


km s"' 






0.599 


0.106 


13.9 


2/orb 


1.197 


0.040 


4.5 


/l 


21.11 


0.040 


4.5 


fl 


30.38 


0.038 


4.0 



6.2. Mode identification witli the FPF method 

We carried out a mode identification of the two detected frequen- 
cies with the FPF method. This method makes use of the fact 
that stellar oscillations cause intensity variations in a rotationally 
broadened line profile due to the Doppler effect. It is capable to 
determine the spherical degree, £, and the azimuthal order, m, of 
a pulsation mode from a time series of line profile variations. 
The contribution of each pulsation frequency to the observed 
line profile variations can be de composed into its Fourier pa- 
rameters, the zero-point, amplitude, phase across the line profile. 
These Fourier parameters are determined from a multi-periodic 
least-squares fit across the line profile, i.e., for each wavelength 
bin, and have a characteristic distribution according to the re- 
lated i and ;7z-values. By comparing the observed Fourier param- 
eters and their counterparts computed from synthetic line profile 
variations, we can determine the pulsational characteristics of 
a pulsation mode. The reduced chi-square values xl of the dif- 
ferent solutions are computed from the complex amplitudes and 
their uncertainties across the line profile. For details we refer to 
IZima, 20061 




10 20 ^ 30 40 

Frequency [d '] 



Fig. 6. Amplitude spectrum of the CoG time-series of the pri- 
mary component. No significant frequency peaks could be de- 
tected in these data. 



The synthetic line profiles were computed assuming that the 
pulsations can be described by a superposition of spherical har- 
monics. Furthermore, we assume a Gaussian intrinsic line pro- 
file, one infinitesimal thin line forming region, and a symmetry 
axis of the pulsation field that coincides with the rotation axis. 

We searched the following parameter space [min; max; 
step]: ( e [0; 15; 1], m e for each (, pulsation ve- 

locity amplitude a e [0; 60; 0.5] km s"' , vsin/e [30;80;1], 
and width of Gaussian intrinsic profile cr E [1; 15; 1] km s"' . 
The inclination value of the RS Cha system is well known 
( [Clausen & Nordstrom, 1980) 1. Assuming that the inclination of 
the rotation axis of the components is identical to that of the sys- 
tem, we set it as fixed parameter to / = 83.4° during the search 
for the optimum in the parameter space. 

Our analysis with the FPF method revealed that both de- 
tected frequencies are high-degree pro-grade pulsation modes. 
Due to the uncertainties of the derived amplitude and phase 
across the line profile, we were not able to unambiguously iden- 
tify a and m but are left with a number of solutions, which still 
will constrain possible theoretical seismic models. The five best 
solutions of each frequency are listed in Table|6l Figure[8]shows 
the Fourier parameters of the observed and best fitting model to 
the two detected frequencies. 

We identified both frequencies as pro-grade modes of high 
degree. The best solution for /i identifies it as sectoral mode of 
t — m - \\, whereas the best identification for fo is a tesseral 
mode with ^' = 10 and m - 6. The uncertainties of this identifi- 
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Fig. 7. Amplitude spectrum of the pixel-by-pixel intensity time- 
series of the primary component. The different panels show the 
data after subsequent pre-whitening of the indicated frequency 
peaks. 



cation are larger for /2 than for /i and in the range of + 1 for both 
quantum numbers. 

We derived a projected rotational velocity of v sin/ =68 + 
2 km s"' which agrees quite well with the previously determined 
value by Alecian et al. (2005) and implies a rotation frequency 
of 0.63 +0.01 d ' . 



6.3. Mode identification w/fh fhe Fourier 2D metinod 

A complementary, and very direct method of analysing the non 
radial pulsation modes present in the LSD-spectra time series 
consists of applying the Fourier 2D method ( |Kennelly, 1994[ 
Kennelly et al., 1996 1, resulting in a representation of the data 
in the frequency - apparent m space. The apparent m is related 
to the structure of the modes present at the stellar surface, with- 
out being identical to the usual azimuthal order m. However, the 
original work by Kennelly, 1994 showed that apparent \m\ scales 



as ( +2 for values close to zero, as { +1 for values lower than 
10 and as i for values above 10. Basically, this method consists 
of doing subsequently a spectral and a temporal Fourier analy- 
sis. All LSD profiles of the time series were normalized to the 
same equivalent width, in order to eliminate the effects of large 
scale variations not related to the signature of pulsations. The 
time average of the mean profiles was then constructed and sub- 
tracted from each mean profile, yielding a time series of normal- 
ized residual profiles. Before doing the 2 dimensional Fourier 
analysis, each residual profile was interpolated on a grid repre- 



Table 6. Mode parameters derived from the appUcation of the 
FPF method. For each pulsation mode, the four best solutions 
are shown, xl is the reduced chi-square value, a is the velocity 
amplitude in km s"' ; cr is the width of the intrinsic Gaussian pro- 
file, both expressed in km s ' . The inclination value was fixed 
at 83.4°; Primary component of RS Cha. 
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senting stellar longitudes, transforming velocities across the line 
profile into longitudes on the stellar equator using the relation 
Av - V sin ; sin <p, where Av is the velocity position within the 
LSD profile with respect to the rest wavelength of the star and 
(p the stellar longitude angle of the star in spherical coordinates 
(under assumption that sin ; is close to unity, which is the case 
for RS Cha). The Fourier transform being very sensitive to gaps 
in the temporal sampling, we selected only continuous series of 
spectra within a night, and interpolated them on a 2" regularly 
spaced temporal grid, where n is selected in order to avoid over- 
or under-sampling. As a result, we obtained F2D spectra for each 
suitable night. In order to take into account the individual qual- 
ity of each night, we decided to combine the spectra of the three 
nights with a time-basis of more than 250 continuously observed 
minutes in the following way: 



F2D 



where aj is a normalized weighting coefficient: 

n ^ 



(1) 



(2) 



k=l 



and S IN j the average signal to noise ratio per spectrum for a 
given night, «spec the number of spectra per night and n the num- 
ber of exploitable nights (in the sense of the Fourier 2D analy- 
sis). The coefficient a^ reflects therefore the total cumulated stel- 
lar photon flux per night, normalized to the cumulated photon 
flux of the full observing run. While this expression is not math- 
ematically rigorous it reflects to some extent the quality of the 
individual nights. 

Figure |9] shows the nightly F2D analysis for the nights of 
12*, 16* and 22"'i of January 2006. [Kurtz & Muller, 1999] es- 
timate the uncertainty in the frequency to be approximately 
l/(4Ar), where AT" is the time span of the data set. In our case for 
a typical observing night of 300 min continuous observations. 



T. Bohm et al.: Discovery of non-radial pulsations in the spectroscopic binary Herbig Ae star RS Cha 







/i 


= 21.11 








0.992 
0.99 






























0.988 
















0.986 














E 


0.984 
































0.982 






****** 









fj = 30.38 d- 



-40 -20 



20 40 



0.996 
0.994 
0.992 
0.99 



0.984 
0.982 
0.98 



■60 -40 



20 40 



0.0007 
0.0006 
0.0005 
0.0004 
0.0003 
0.0002 
0.0001 



■ TItT 



0.0008 
0.0007 
0.0006 
^ 0.0005 
I 0.0004 
I" 0.0003 
" 0.0002 - 
0.0001 

- 
-0.0001 









data 1 — 1 — 1 








1=11, m=11 

1=11, m=9 




a 
























■-. 


































-40 -20 20 40 

velocity [km s' j 











data 












kOi IT 


=5 














- 

'Ml] 























































-60 -40 



20 40 



Fig. 8. Fourier parameters across the line profile of the detected pulsation frequencies. Each panel shows from top to bottom: zero- 
point, amplitude and phase. For both frequencies, the two best fits to the observed zero point, amplitude in units of the continuum 
and the phase distribution in units of 2n are shown. The uncertainties are the formal errors computed from the least-squares fitting. 
Primary component of RS Cha. 



Table 7. Frequency search and mode identification. Summary of 
the weighted F2D analysis on the primary component. Columns 
are as following: (1) frequency (error bars are of the order of 
1.2d"') (2) possible identification with FPF result, (3) apparent 
\m\, (4) expected i value. 



d-l 


FPF frequency app. \m\ 


e 


20. 


/, 8-10 


8-10 


24. 


4 


2 or 3 


12. 


8 


8 


4. 


2 






this would yield an uncertainty of 1.2d in the frequency deter- 
mination. This, together with large-scale beating effects on the 
pulsational signature and significant S/N ratio variations from 
one night to another might account for slightly different po- 
sitions of frequencies in the F2D space. However, the "gen- 
eral" pattern is preserved from night to night, as can be seen 
in Figure |9l and the weighted F2D periodogram in Figure [10] 
exhibits the strongest recurring oscillation modes. 

Table |7] summarizes the results for the primary component 
based on the F2D analysis. 

The F2D method confirms the existence of the high degree 
mode /i determined by the FPF method (/i = 21.11 d"') and 
indicates a degree { of 8-10, within the error bars of the FPF de- 
termination. While single night F2D periodograms show some 
high frequency signal around 29. d ' (e.g. night of Jan 12*), 
which could correspond to /2 (FPF), this is not confirmed sys- 
tematically. A low frequency signal around 4. d"' is very likely 
a fundamental oscillation frequency due to the window function 
(tohs typically 1/4 d) and is present in all F2D periodograms seen 
in Figure |9l Other frequencies found in the weighted F2D peri- 



odogram around 12. d ' and 24. d ' (with different apparent \m\ 
values) have not been revealed by the FPF method. 

7. Pulsational analysis of the secondary component 

The analysis of the secondary component was performed on a 
selected data subset, following the same selection rules as for the 
primary component. This procedure left 165 spectra that were 
suited for the further analysis. Also for the second component, 
we had to cut out the central parts of the line profile between 
+70 km s"^ in Doppler velocity. Therefore, we also had to use 
the CoG of the profile to estimate the radial velocity variations. 
We carried out the frequency search and mode identification in 
the same manner as for the primary component. 

7. 1 . Frequency analysis of the secondary component 

The Fourier spectrum of the CoG time series (see Figure [TTI ) 
shows a dominant frequency peak with an amplitude of 1 km s ' 
at twice the orbital frequency 1.19 d ' . After pre- whitening of 
this frequency, the low frequency domain below 6 d ' does not 
show any more significant frequencies and has a much lower 
noise level than the primary component. Strong 1 d ' -aliasing 
complicates the detection of further frequencies. The highest re- 
maining peak is at 12.8 d"' with an amplitude of 0.78 km s"' , 
but the one-day aliasing peak at 13.9 d"' has a similar ampli- 
tude. A similar situation occurs at 20.1 d"' . We selected two 
peaks that had a S/N above 4 and that resulted in the lowest 
residuals after pre-whitening, /i = 12.81 and /2a = 20.11 d"' . 
After pre-whitening with these frequencies, no further signifi- 
cant peaks can be found in the Fourier spectrum. 

The analysis of the PIV across the line profile (see Figure[T2]i 
revealed a dominant peak at 0.19 d ' , which is a one-day alias- 
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Fig. 10. Weighted F2D periodogram of the Primary component 
of RS Cha for the nights of Jan 12*, 16* and 22"'' 2006. 



apparent rr^ 

2D Fourier analysis RS Cha Primary 22 jan 2006 




Fig. 9. F2D analysis on the Primary component of RS Cha for 
the nights of Jan 12*, 16* and 22"'' 2006. Apparent |m| scales as 
{ +2 for values close to zero, as ^ +1 for values lower than 10 
and as £ for higher values ( Kennelly, 1994[ i. 



ing of twice the orbital frequency at 1.19 d"' . Also the orbital 
frequency had a significant amplitude in the PIV which, as has 
been mentioned before, is a result of imperfect orbital fitting. 



Table 8. Frequencies detected in the pixel-by-pixel intensity 
variations of the secondary component of RS Cha. The value 
Apiv is the integral of the pixel-amplitude across the line pro- 
file between -70 and 70 km s"' . For each frequency, the S/N is 
computed for the wavelength bin where its amplitude is highest. 



ID 


freq. 


Apiv 


S/N 




d-' 


km s"' 




/orb 


0.599 


0.060 


5.7 


2/orb 


1.197 


0.138 


8.2 


/l 


12.81 


0.066 


3.9' 


f2b 


19.11 


0.061 


5." 


h 


24.56 


0.054 


4.0 



We pre-whitened the data with these two frequencies and found 
fx - 12.81 d"' as highest remaining peak. After further pre- 
whitening, the highest peak is at f2b - 19.11 d"' which is a 
one-day aliasing of f2a- Pre-whitening with fy, resulted in much 
lower residuals than pre-whitening with /2a, which contradicts 
the frequency analysis of the CoG. Since we used the PIV for 
the mode identification with the FPF method, we accepted /a;, as 
the 'true' frequency. After pre-whitening, one more frequency at 
/3 - 24.56 d"' with a S/N of 4 was detected. No other significant 
peaks could be extracted from the data. The detected frequencies 
are Usted in Table|7.1] 



7.2. Mode identification witli fhe FPF metliod 

We carried out a mode identification with the FPF method for 
the frequencies /i, f2b, and fj, and searched the same parame- 
ter space as for the primary component. The results are listed in 
Tabled According to our mode identification, fi is a low-degree 
pro-grade mode, whereas the two other higher frequencies are 
high-degree pro-grade oscillation modes. In comparison to the 
mode identification of the primary, the values of the goodness of 
fit (xl) are higher for f2b and fi , which implies a lower reliability 
of the derived identification. For these two frequencies, the de- 
gree, (, is much better constrained than the value of m. Figure[T3] 
shows the Fourier parameters of the observations and the two 
best fitting models for each mode. 
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Fig. 13. Fourier parameters across the line profile of the detected pulsation frequencies. Each panel shows from top to bottom: zero- 
point, amplitude and phase. For both frequencies, the two best fits to the observed zero-point, amplitude in units of the continuum 
and the phase distribution in units of Itt are shown. The uncertainties are the formal errors computed from the least-squares fitting. 
Secondary component of RS Cha. 



Table 9. Mode parameters derived from the application of the 
FPF method. For each pulsation mode, the four best solutions 
are shown, xl is the reduced chi-square value, a is the velocity 
amplitude in km s"^ ; vsin/ is the projected rotational velocity; 
cr is the width of the intrinsic Gaussian profile, both expressed 
in km s"' . The inclination value was fixed at 83.4°. Secondary 
component of RS Cha. 
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Table 10. Frequency search and mode identification. Summary 
of the weighted F2D analysis on the primary component. 
Columns are as following: (1) frequency, (2) possible identifica- 
tion with FPF result, (3) apparent \m\, (4) expected { value. 



d-' 


FPF frequency 
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19. 


f2l, 
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12. 
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7.2.1 . Results with the Fourier Fourier 2D method 

Figure [14] shows the nightly F2D analysis for the nights of 12*, 
16* and 22"'' of January 2006. The weighted F2D periodogram 
in Figure [T5]clearly exhibits the presence of 2 main frequencies 
around 12 d ' and 19 d ', which confirms the frequency deter- 
mination by the FPF method for f\ and /a;,. The F2D method 
indicates an ( value of 2 or 3 for f\ and confirms the FPF de- 
termination (within the error bars). However, the F2D method 
attributes a low I value of to f2b and identifies therefore this 
mode as radial, quite in contradiction with the FPF result for the 
same frequency, indicating a high degree ( = 9-13. 

Table [TO] summarizes the results for the primary component 
based on the F2D analysis. 



8. Discussion and conclusions 

For the first time we detected by direct spectroscopic means non- 
radial pulsations in both components of a binary Herbig Ae star 
We identified with a high level of confidence in the profile varia- 
tions of the primary and secondary component two, respectively 
three pulsation frequencies. A first identification by two comple- 
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Fig. 11. Amplitude spectrum of the CoG time-series of the sec- 
ondary component. The different panels show the data after sub- 
sequent pre-whitening of the indicated frequency peaks. 



Fig. 12. Amplitude spectrum of the pixel-by-pixel intensity time- 
series of the secondary component. The different panels show 
the data after subsequent pre-whitening of the indicated fre- 
quency peaks. 



mentary methods revealed very different main pulsation modes 
for the two components: while the dominant mode of the pri- 
mary component seems to be a high degree pro-grade mode with 
^ = 10 or 11, the dominant mode of the secondary component is 
identified as a low degree mode with ^ = 0,1 or 2. Two differ- 
ent methods have been applied to perform a mode identification 
of both components. The Fourier 2D method as presented by 



Kennelly, 1994 provides direct insights on the spatial structure 
of the oscillation. However, the derived apparent \m\ value seems 
to be linked to { in the way described here-above, but more de- 
tailed simulations for different configurations have to be carried 
out in order to understand potential biases in the degree I deter- 
mination. We therefore privilege the E determination as provided 
by the FPF method. TablefTTIsumm arizes the main results for the 
dominant pulsation modes of each component. As can be seen, 
both methods yield very different results on the determination of 
the degree t of /a;, of the secondary component, perhaps due to 
an insensibility of the F2D method to tesseral mode signatures. 

The precise redetermination of the binary orbit reveals a 
slight gradient change of the observed - calculated phase shift 
curve, which might be the signature of changes in the orbital pe- 
riod of the order of AP/P = (7.5 + 0.7) lO"*". This tendency, if 
confirmed, might be due to the presence of a third body in the 
RS Cha system or due to intrinsic orbital period changes, which 
then need to be explained. 

Next steps require a precise redetermination of the funda- 
mental parameters and a search for the intrinsic rotation periods 
of both components. The first frequencies are now known and 



their associated modes are constrained by the present work. In 
subsequent work, a numerical model will be developed in order 
to constrain the internal structure of both components. 



Table 11. Frequency search and mode identification. Summary 
of the main results for the dominant frequencies of each 
component. The {E,m) couples indicated for the FPF method 
correspond to the best identification, in decreasing order For 
each of the frequencies detected by the FPF method, we identi- 
fied an I value as detected by the F2D method. 
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Fig. 14. F2D analysis on the Secondary component of RS Cha 
for the nights of Jan 12*, 16* and 22"<^ 2006. 
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Fig. 15. Weighted F2D periodogram of the Secondary compo- 
nent of RS Cha for the nights of Jan 12*, 16* and 22"'' 2006. 
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